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Recent studies have shown that the primordial non-Gaussianity affects clustering of dark matter 
halos through a scale-dependent bias and various constraints on the non-Gaussianity through this 
>i-s ' scale-dependent bias have been placed. Here we introduce the cross-correlation between the CMB 

I , lensing potential and the galaxy angular distribution to effectively extract information about the 

f~>) ■ bias from the galaxy distribution. Then, we estimate the error of non-linear parameter, /nl, for the 

^vj ' on-going CMB experiments and galaxy surveys, such as Planck and Hyper Suprime-Cam (HSC). 

We found that for the constraint on /nl with Planck and HSC, the wide field galaxy survey is 
preferable to the deep one, and the expected error on /nl can be as small as: A/nl ~ 20 for feo = 2 
, and A/nl ~ 10 for ho — 4, where bo is the linear bias parameter. It is also found that future wide 

^Lf ' field galaxy survey could achieve A/nl ^ 5 with CMB prior from Planck if one could observe highly 

, -»^ , biased objects at higher redshift (2 ~ 2). 

Q ! I. INTRODUCTION 

u: . . 

The cosmic microwave background (CMB) temperature anisotropy is a quite useful probe for cosmology. The 

contribution to the anisotropy is dominated by fluctuations at the last scattering surface. The CMB photon, however, 

1 '• encounters the large-scale structure along the line of sight and some additional effects are imprinted on the temperature 

O ' and polarization as secondary anisotropics. The deflection of the CMB photon due to gravitational potential produced 

-*— > , by the large-scale structure is one of them. The effect of the gravitational lensing on the CMB photon through the 

^ ' large-scale structure is known as the CMB lensing [1| . An on-going CMB observation by Planck [2| or various ground- 

I— I , based experiments are expected to detect this signal, while the effect of the CMB lensing are imprinted on small scales 

which the Wilkinson Microwave Anisotropy Probe (WMAP) satellite could not resolve. The lensing effect reflects the 

vj ] late time evolution of the universe at relatively low-redshifts. Therefore, the lensing information plays an important 

Jr> i role to determine the cosmological parameters, such as the neutrino mass, the cosmological constant, the equation of 

Qs ' state parameter of dark energy and so on. 

^+ . The large-scale structures arc formed at relatively late time and they become the source of the gravitational 

C^ ' potential. They are correlated with the CMB temperature anisotropy through the Integrated Sachs- Wolfe (ISW) 
effect, which generates the secondary anisotropics due to the time variation of the potential Q. The cross-correlation 
^S ' has an advantage for observations of ISW effect whose signal is weak. Cross-correlations with complementary probes 
f^ [ are expected to provide additional information on top of their respective auto-correlations. Similarly, we expect that 
lensing potential should correlate with the large-scale structure and the information from their cross-correlation may 
precisely determine the cosmological parameters. 






X 



Recently, the deviations from Gaussian initial conditions (primordial non-Gaussianity) arc intensively focused on 
and discussed. Inspection of them offers an important window into the very early Universe because non-standard 
models of inflation allow for a large non-Gaussianity while standard single-flcld slow-roll models predict the small 
deviations from Gaussianity. The most popular method to detect the primordial non-Gaussianity is to measure 
higher-order correlations of CMB anisotropics and distributions of galaxies, for example, the bispectrum or the three- 
point correlation function of CMB 0, Q or the large-scale structure bispectrum @-Q . 

Some studies have shown that the primordial non-Gaussianity affects cluste ring of dark matter halos through a scale- 
dependent bias, both by analytic calculations and by N-body simulations |10l - [l4| . By considering the scalc-dcpcndcnt 
bias, one can constrain on the non-Gaussianity by the powe r sp ectrum. Various constraints on the non-Gaussianity 
through the scale-dependent bias have already been placed [14| - [l9| . One may expect, however, that a constraint on 
non-linear parameter /nl, which describes the primordial non-Gaussianity, is degenerated with other cosmological 
parameters and has a large error. Because of such degeneracy, especially with the linear bias bo, it is important to 
combine unbiased observations which are sensitive to the matter power spectrum in the near universe, like CMB 
lensing, shear and so on. Here we use the cross-correlation between the CMB lensing potential and the large-scale 
galaxy distribution and estimate expected errors of non-linear parameter /nl for future CMB experiments and galaxy 
surveys. We expect that this cross-correlation may play a role to break degeneracy and give us more stringent 



constraints. 

This paper is organized as follows. We review the scale-dependent bias due to the non-Gaussianity in section [II] 
and the theory behind the cross-correlation between CMB lensing potential and galaxy distribution in section IIIII 
In section IIVI we describe the survey model for the galaxy distribution in the Hyper Suprime-Cam (HSC) survey, 
which is a fully funded imaging survey at Subaru telescope. In section IVl and section IVTl we explain the method of 
our analysis. Finally, in section fVIII and I Villi we discuss the results and summarize our conclusions. Throughout this 
paper wc assume a spatially flat universe for simplicity. 

II. SCALE-DEPENDENT BIAS 

Deviations from Gaussian initial conditions are commonly parameterized in term of the dimcnsionlcss /nl parameter 
and primordial non-Gaussianity of the local-type is defined as |J| 

$ = + /nl(<^'-(0')), (1) 

where $ denotes Bardeen's gauge-invariant potential and (j) denotes a Gaussian random field. On subhorizon scale, 
$ = — ^, where '^ denotes the usual Newtonian gravitational potential related to density fluctuations via Poisson's 
equation. For example, simple slow- role inflation gives a parameter /nl of the order of 10^^ — 1 , 20i-,22| . On the other 
hand, large values of /nl can be expected in models of multifield inflation, tachyonic preheating in hybrid inflation 
j23l | or ghost inflation [24|, for instance. Thus, the information about the inflation physics is closely related to the 
parameter /nl- 

Recent studies show that the effect of the primordial non-Gaussianity of the local-type is seen in the clustering of 
halos through a scale-dependent bias, 

Pg{k) = blP{k) ^ [bo + A6(fc)]2p(fc), (2) 

where Pg{k) and P{k) are the power spectrum of galaxy and matter density fluctuations as a function of the wave 
number k, respectively. Bq is the Gaussian-case bias which relates the galaxy density fluctuations with the matter 
density fluctuations, and Ab{k) represents the scale-dependence due to the non-Gaussianity [IOtIJIi 

^^^''> - D{z)k^T{k) ' ^'^> 

where D{z) and T{k) are the growth rate and the transfer function for linear matter density fluctuations, respectively. 
5c — 1.68 is the threshold linear density contrast for a spherical collapse of an overdensity region. Primordial non- 
Gaussianity of the local-type gives rise to a strong scale-dependent bias on large scales, while the bias is roughly 
constant on large scales in the Gaussian case. However it is necessary to emphasize that the constraint through 
the scale-dependent bias is sensitive only to the local-type non-Gaussianity. For the constraints on the other non- 
Gaussianity models we must consider higher-order correlation such as bispectrum, trispectrum and so on. 

III. THE ANGULAR POWER SPECTRUM 

The cross-correlations, for example, between CMB and galaxy, are well known as providing additional information 
other than their respective auto-correlation. In Ref. [25[, they investigated the cross-correlation between the shear 
of CMB lensing and halos. In this paper, we introduce the cross-correlation between the CMB lensing and galaxy 
angular distribution to estimate errors in constraining cosmological parameters. 

A. Galaxy Distribution 

Probably the most obvious tracers of the large-scale density field in the linear regime are luminous sources such as 
galaxies at optical wavelengths and AGNs at x-rays and/or radio wavelengths. The projected density contrast of the 
tracers can be written as 

5g{n) = / dz—6g{xn,z), (4) 



where 5g represents the density contrast of tracers, n is the direction to the hne of sight, dN/dz is a normahzed 
distribution function of tracers in redshift such that J dzdN/dz = 1 and xi^) is the comoving distance to the redshift 
z. We assume the following analytic form of the normalized galaxy distribution function, 



dN 
dz 
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exp 



(5) 



where a, /3 and zq are the free parameters. In this parameterization a and (3 denote the slope of the distribution 
at low and high-redshifts, respectively, and zq determines the peak of the distribution. We assume that the tracer 
density field is related to the underlying matter density field via a scale- and redshift-dependent bias factor, so that 
5g{k,,z) = h{k,z)5{k,,z). On large scales, where the mass fluctuations are small (5 <^ 1, the perturbations grow 
according to the linear growth rate, (5(fc, z) = d{k)T{k)D{z), where 6{k) is the primordial value of matter density and 
T(fc) is the transfer function. The linear angular power spectrum of the galaxy distribution for a flat universe is given 
by 



cr 



^ ' edkp{k)i\f{kf 



(6) 



where 



Af (fc) = J dz^b{k, z)T{k)D{z)ji{kx). 



(7) 



and P{k) is the linear power spectrum as a function of the wave number k and ji{kx) is a spherical Bessel function. 
In order to estimate errors in parameters and signal-to-noise ratios, we need to describe the noise contribution due 
to the finiteness in numbers of sources associated with source catalogs. We can write the shot noise contribution as 






(8) 



where til is the surface density of sources per steradian and related to the total number of available samples iVg as 
riL = iVg/47r/sky. We show the angular power spectrum of the galaxy distribution, Cf^, and the noise spectrum, N^^, 
in the left panel of Fig. [TJ 



B. CMB Lensing Potential 



We consider the potential that dcfiects CMB photons. The relationship between the lensed temperature anisotropy, 
T{n), and unlcnscd one, T{h), is related by T{h) = T{h + d) and the deflection angle d[n) is related to the line of 
sight projection of the gravitational potential ^(xw,?]) as d{fi) = Vip{h), where 



V'(n) = -2 



lK{X*)fK{x) 



x). 



(9) 



Here '<p{n) is the lensing potential, JkIx) is the angular diameter distance, x is the radial comoving distance along the 
line of sight, and x* denotes the distance to the last scattering surface. For a flat universe angular diameter distance 
is related to the comoving distance as /i<-(x) = X- The angular power spectrum of the lensing potential for a fiat 
universe can be written as 
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(10) 



where 



Af(fc) 



dxT^ik;rio - x) 



x*x 



jiikx)- 



(11) 



In the linear theory, we define a transfer function for the gravitational potential T<i,{k;'q) so that P'i,{k\ri) = 
Ti{k;r^)P{k). 

The lensing potential can be reconstructed using quadratic statistics in the temperature and polarization data that 
are optinfized to extract the lensing signal. To reconstruct the lensing potential f/;, one needs to use the non-Gaussian 



information imprinted into the CMB. Lensing conserves surface brightness, so that the probabihty distribution function 
of the temperatures remains unchanged. Therefore the lowest order nonzero estimator of the lensing potential is 
quadratic. This quadratic estimator has been investigated by [2fl, 123] and the minimum variance estimator was given 
by [28| . A quadratic estimator in the flat-sky approximation generally has the form [2g| 

m) = NiL) I ^em'iL - l)g{l, L), (12) 

where and G) are lensed temperature and/or polarization modes on the sky, i.e., Q, Q ~ T, E, B. The optimal 
weight g(l,L) and normalization N{L) for each mode are found using the fact that the deflection position can 
be written as a first order expansion of the displacement around the undeflected position, 9(n) = &{ri + d) = 
Q(n) + V*^'(n)Vi0(n). Requiring the estimator to be unbiased and minimizing the variance, the optimal weight for 
TT estimator is 

^ {L~l)-Lqr.-ii+l-LQ 

^^ ' ' or'tot/Vtot ' ^ ' 

'^^l ^\L-l\ 

where C; (C/) is the unlensed (lensed) temperature power spectrum. For other estimators, C; (C;) represents the 
temperature or polarization one. The superscript "tot" originates from the fact that the lensed CMB and the noise 
enter in the variance, (7*°* = Ci + Ni. 

With the definition in Eq. (1121) . the lowest order noise of the lensing reconstruction equals to the normalization 
which is determined by 



SiO)mL)\')^NiL) 



J (0^ [^^ " '^ ■ ^^l^-'l + ' ■ ^^'] "" ^^'' ^^ 



-1 



(14) 



(2^) 
Physically the variance is a combination of the noise introduced by primary anisotropics themselves and the instru- 



mental noise. The all-sky generalization is presented in Ref. [27 



Here, the noise power spectrum of the CMB experiment reads 

Ni"^^^ = iOpwHMAxf exp [l{l + 1)0|whm/81ii2] , (15) 

with X G {T, E, B}, where Ax is the temperature and polarization sensitivities per pixel of the combined detectors 
and 0FWHM describes the spatial resolution of the beam. These values are given for each frequency bands v and we 
show the values for some CMB experiments in Table ID When there are multiple frequency bands or channels, the 
global noise of the experiment is given by 



Arf^ = 






(N^'^ 



(16) 



where the sum is over the individual channels. We show the angular pawer spectrum of the CMB lensing potential, 
C[ , and its noise spectrum, A^^ , for various CMB experiments in the right panel of Fig. [TJ As Planck does not have 
much sensitivity to reconstruct the lensing potential from the polarization components, TT provides the best estimator 
for the Planck. For the reference experiment like the CMBPol, the lensing potential, however, is reconstructed from 
polarization components and EB provides the best estimator. 

C. Cross-Correlation: Galaxy &: Lensing Potential 

We focus on the linear cross-spectrum of the galaxy with the CMB lensing potential, 

Cfs = - f fc2dfcp(fc)Af (fc)Af (fc). (17) 

The most important assumption we have made so far is that the galaxy distribution and the lensing potential is linear 
and Gaussian. On small scales this will not be quite correct due to non-linear evolution. For simple models, fits to 



numerical simulation like the HALOFIT code of Ref. |3l| can be used to compute an approximate non-linear power 



Experiments 


/sky 


V [GHz] 


6'fwhm 


At 


Ap 


Planck [2] 


0.65 


100 


9.5 


6.8 


10.9 






143 


7.l' 


6.0 


11.4 






217 


5.0 


13.1 


26.7 


PolarBear [29] 


0.03 


90 


6.7 


1.13 


1.6 






150 


4.0' 


1.70 


2.4 






220 


2.7 


8.00 


11.3 


CMBPol [30] 


0.65 


100 


4.2 


0.87 


1.18 






150 


2.8 


1.26 


1.76 






220 


1.9' 


1.84 


2.60 



TABLE I: The current designs of CMB experiments, ^fwhm is the Gaussian beam widtli at FWHM, At and Ap are the 
temperature and polarization noises, respectively. Planck and CMBPol are the satellite experiments and PolarBear is the 
ground based experiment. 
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FIG. 1: (Left) Angular power spectrum of the galaxy distribution, Cf®, for the Gaussian initial condition (/nl = 0), and the 
galaxy shot noise, l/riL. We show Cf® for various forms of galaxy sampling model [Eq. J5|], (a, /?) = (0.5, 3.0) (solid line) or 
(a, /3) = (2.0, 1.5) (dashed line), and zq = 0.8 (thin line) or 2o = 1.2 (thick line). The dotted, dot-dashed and dash-dot-dashed 
lines show shot the noise contributions from til = 100, 200, 400 [deg~^], respectively. (Right) Angular power spectrum of CMB 
lensing potential, Cf (solid line), and lensing reconstruction noise, N"^^ (non-solid lines). Each of the noise power spectrum 
indicates those for Planck (dashed line), PolarBear (dot-dashed) and CMBPol (dot-dash-dotted line), respectively. 

spectrum. A good approximation is simply to scale the transfer functions T(/s) of Eq. (O, pT|) so that the power 
spectrum has the correction from the non-linear effect 



T(fc) 



T{k). 



Ipn 
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(fc) 



p(fc) 



(18) 



We also include the other cross-correlation components, Cf^, C^ and Cj ®, for the estimation of the parameter 
errors. However, we assume that there is no cross-correlation between the polarization and the lensing potential or 



the galaxy distribution, C, 



£■!/> 



cf' 



0. This is because the polarization is mainly produced by the Thomson 



scattering at the last scattering surface while the lensing potential and the galaxy distribution exist in the late-time 
universe. We show the angular power spectrum C^ in Fig. [51 The redshift dependence and the effect of the primordial 
non-Gaussianity through a scale-dependent bias are clearly seen in that figure. 
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FIG. 2: The cross-correlation power spectrum between the CMB lensing potential and the galaxy distribution, l{l + 1)C; /27r 
for /sky ~ 0.1 and A'^g — l(f . The solid, dashed, dotted and dot-dashed lines correspond to /nl ~ 0, ±50, ±100, ±500, 
respectively. 



IV. MODELING GALAXY SAMPLE 



We showed the analytic form of the normalized galaxy distribution function in Eq. ([5|). The mean redshift Zm is 
related to the peak redshift zq a-nd determined by 



= / dzz 



dN _ zgT [{a + 2)//3] 

dz r[(a + i)//?] ■ 



(19) 



Jab limit 



i{S/N = 10) 



g{S/N = 5) 



r{S/N = 5) 



z{S/N = 5) 



22.97 
23.84 
24.22 
24.81 
25.04 



1 mins. 

5 mins. 
10 mins. 
30 mins. 
45 mins. 



3 mins. 
15 mins. 
30 mins. 
90 mins. 
130 mins. 



1.1 mins. 
7 mins. 
12 mins. 
34 mins. 
50 mins. 



0.3 mins. 

1.4 mins. 

3.5 mins. 
8.1 mins. 
13 mins. 



TABLE II: Exposure time for the bands, i,q,r,z. The relation between magnitude limit and exposure time was scaled for the 
published Subaru Suprime-Cam specification [33. [35l|. 

The relation between zq and Zm is, for example, zq = Zm/0.64 for (a, (3) = (0.5, 3.0) and zq = Zm/1.41 for (a, (3) = 
(2.0, 1.5). In this paper, we consider a wide field survey such as the on-going Hyper Suprime-Cani (HSC) project. 
This is a fully funded imaging survey at Subaru telescope. The surface density ^l and the mean rcdshift z^ are 
related to the exposure time icxp as [321 133| , 



0.9 



/ tc 



\ 30min 



0.067 



0.44 



riL = 35 — ™P_ j [arcmin 

\30min/ 



(20) 
(21) 



In Ref. [S^l and [33[, (a, /3) = (0.5, 3.0) and (2.0, 1.5) are adopted, respectively. In this paper, we adopt both cases 
and compare the differences between the survey models. 

The validity of the above form of the galaxy distribution is shown in Ref. [33|. They compared it with the 
Canada- France-Hawaii telescope (CFHT) photometric redshift data [SJl- The relationship between magnitude limit 
and exposure time was scaled for the published Subaru Suprime-Cam specification |35j . and these data are shown in 
Table HIl for the i,g,r,z passbands. 

The total survey area can be expressed as [33 



/ field of view 






Ttotal 



1.1 X tcxp + to 



(22) 



where we assume that the field of view is 1.5°, the total observation time Ttotai is fixed as 800 hours, and the overhead 
time is modeled by constant, top = 5 min, plus a fraction (10%) of the exposure time texp for one field of view. 

V. THE FISHER MATRIX ANALYSIS 

For our Fisher matrix analysis, we refer to the method of Ref. [36| and expand it to take into account the cross- 
correlation between the lensing potential and the galaxy distribution. In Ref. [3Q|, the 5x5 covariance matrix is 
calculated for primary CMB and CMB lensing. In our case, we expand it into 8x8 covariance matrix for the 
cross-correlation between CMB lensing and galaxy distribution. 



Likelihood Function 



Each data points have contributions from both signal and noise. If we assume both contributions arc Gaussian 
distributed, we can write the likelihood function of the data given the theoretical model as 



'^<''i'"«;^^"-(-5'"P<«)-i<')' 



(23) 



where d ~ i'^im^^im^'^inn^hn) ^^ ^^^ ^^^^a vector, Q — (6*1,^2, ■•■) is a vector describing the theoretical model 
parameters, and C{Q) is the theoretical data covariance matrix represented by both signal and noise. For it to be a 
good estimate, we would like it to be unbiased, i.e., (0) = Gq , where Qq indicates the true parameter vector of the 
underlying cosmological model, Q is the one constructed by the data vector d they minimizing the likelihood function 
C{d\<d) {i.e., the so-called best fit model) and (. . .) denotes an average over many independent realizations. 



We can derive the effective chi-square, xls = Sxy Sim (^2) ln£, from (|23|) as 

Xe\-E(2' + l)(p+l"l'^l)' (24) 

where J2xy represents the summation for each modes, X,Y = {T,E,(t),g}, and \C\ denotes the determinant of the 
theoretical data covariance matrix, 

\C\ = {Cf''f{Cf'f-Crcr{Ct'f-Cf^Cf{Cf''f 



(26) 



Here D is defined as 



D = +2(7P(7r^(Cf^)2 + 2((7r^)2cf^q^9 



_^TT^EEf^4>g\2 _ ^TT/^EE/^4>g\2 __ nr"^"^ C^^ C^^ C^^ 
_^'ipiiQggf^TE\2 _ Qip^ ^gg r^TE\2 _ 2(j^^ ^99 (jTE ^TE 
_^EE^-4'iilQTg-.2 __ QEE^ii-4'/^Tg-.2 _ 2qEE ^^i^i ^Tg ^Tg 
_^EE^gg/^T^-.2 _ ^EE^ggi^Ti'^ _ 2C^^ C^^ C^^ C^^ 
,^TT^EE^ii-4'^gg , ^TT^EE^tptl'^gg , ^TT^EE^i'tl' ^gg , ^TT^EE^i'ii^gg 

L L t I, L L i t, L L t t L L l t 

j^2C^^ C^^ C^9 ^^g ^ 2C^^C^'^C^^C'^'^ + 2C^^C^'^C^^C'^^ + 2C^^C'^'^C'^^C''^^ (27) 

i L L L t' L L L f' L L L f- L L L ^ ' 

In the above expression, we have assumed that the polarization component does not correlate with the lensing potential 
and galaxy distribution, so we put C; = C; ^ = 0. 

On the other hand, the mock data covariance matrix C is given from the simulations and defined as C* = {dS). 
We can estimate the power spectrum of the mock data through the following definition, 

^aL*aL = (2/ + l)Q^^. (28) 

From Bayes' theorem, we assume ([24|) to be the distribution of theoretical data covariance matrix (7(6) when mock 
covariance matrix C is given. Then, we can account C to be a variable and C to be a constant. All expressions 
introduced so far assume a full sky coverage survey. However, real experiments can only see a fraction of the sky. We 
introduce a factor /sky, where /sky denotes the observed fraction of the sky in the effective x^- We are interested only 
in the confidence levels, so the normalization factor in front of the likelihood function ([23]) is irrelevant. We normalize 
as x^g = if C* = C* by adding arbitrary constant and redefine x^g from ([24]) as 

Xe'ff = > '(2? + l)/sky|^+ln^-4|, (29) 




and \C\ denotes the determinant of the mock (observed) data covariance matrix, 



B. Fisher Information Matrix 

The Fisher matrix formalism can be used to understand how accurately we can estimate the values of vector of 
parameters for a given model from one or more data sets [33|- The Fisher matrix approximates the curvature of 
the likelihood function L around its maximum in a space spanned by the parameters 9. The usual formula requires 
a slight generalization to account for the possibility that different surveys may only partially overlap in sky coverage 
as we shall show below. The likelihood function should peak at ~ 8o, and can be Taylor expanded to second order 
around this value. The relevant term at second order is the Fisher information matrix, defined as 






Fv ^ -i^TTinr ■ (31) 



eo 



From the Cramer- Rao inequality, the marginalized error on a given parameter 0i is given by cr{9i) = \/{F^^)ii for an 
optimal unbiased estimator such as the maximum likelihood. 

Substituting equations (|23|) and (|29|) into the above expression, the Fisher information matrix is written by 



'max an^x anyy 



1=2 XX'.YY' 



de. 



IXX'YY' 



de. 



(32) 



where i, j run over the cosmological parameters, /„iax is the maximum multipole available given the angular resolution 
of the considered experiment, and XX', YY' e {TT,EE,TE,iljip,Tilj,gg,Tg,ipg}. The matrix Cov; is the power 
spectrum covariance matrix at the /-multipole, 

f ^TTTT ^TTEE ^TTTE ^TTcj)cj) ^TTgg ^TTTc/) ^TTTg ^TTipg \ 

'^ttee '^eeee '^teee 

Ettte '^teee '^tete 
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ExTTg St(,V'^ ^TgTg ^T^-Tg '^Tg^g '^Tg^.g 

\ ^TT'ipg ^^^^g !^ipggg ^Ttp-ipg ^Tg^g ^ipgipg 



Covi 
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where the auto correlation coefficients are given by 



^TTTT 
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(33) 
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(34b) 

(34c) 

(34d) 

(34e) 
(34f) 

(34g) 
(34h) 
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while the cross-correlation ones are 
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VI. FORECASTS 



We estimate the parameter errors for Planck satellite using Fisher analysis following the method introduced in Sec. 
IVl Our fiducial cosmology is based on the WMAP 7- year result |33| within a flat, ACDM model framework. 
The fiducial model parameters we consider arc given as : 



e = { rife/i^, ^ch^, Qa, t, f^, Fho, ns, A^(fco) x 10^, wq, N^s, as, /nl, ^o } 
= {0.02258, 0.1109, 0.734, 0.088, 0.02, 0.24, 0.963, 2.45,-1.0, 3.0, 0, 0, 2.0}, 



(36) 



where fib, ^c and ilA are the density parameters for baryon, cold dark matter and dark energy, respectively, h is 
the Hubble constant, r is the Thomson scattering optical depth to the last scattering surface, fi, is the mass density 
of the massive neutrino relative to the total matter density: f^, = $7^/57,^, Iho is the primordial helium fraction, 
Ug is spectral index of the primordial power spectrum, Ap(fco) is the amplitude of the primordial power spectrum 



TABLE III: Survey parameters for our calculation. 
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normalized at fcg = 0.002 Mpc^^, w is the equation of state parameter of dark energy, N^i^ is the effective number of 
neutrinos, as is the running index, /nl is the non-linear parameter which represents the primordial non-Gaussianity 
and 6o is the linear bias parameter. Because we assume a flat universe the Hubble parameter is adjusted to keep 
our universe flat when we vary ^a- For neutrino parameters, we assume the standard three neutrino species. In 
our analysis, the non-linear parameter /nl and the linear bias parameter 6o are determined by the galaxy surveys 
only, and the CMB experiment plays a role in breaking the parameter degeneracies. We use CAMB code [23 and 
HALOFIT code l3l| to calculate the angular power spectrum C^^ and the non-linear region of the angular power 
spectrum of the galaxy distribution and lensing potential. 

In our estimation, we include the information from temperature anisotropics, E-mode polarization and reconstructed 
lensing potential. The range of multipoles are 2 <l < 2500 for Cf"'^ and Cj^^ and 2 < I < 1000 for Cf, respectively, 
and survey area is taken as f^^ = 0.65 for CMB survey. On the other hand, we include the information from 

galaxy survey, Cf^, where the range of multipoles are 2 < I < 1000 and survey area is f^-^y'^^^ = 0.10. We assume 
that there is no correlation between different patches, so that the area where there is a correlation between CMB 
and galaxy survey corresponds to the galaxy survey area /galaxy We summarize the values we used mainly in the 
following calculation in Table IIIIl 

The non-linear effect on the angular spectrum of the galaxy distribution begins to appear at I > 100. As the 
calculation of Fisher matrix assumes that all fields are random, the non-linear region of the galaxy distribution is 
inadequate for this calculation. However, because the auto correlation signal of the galaxy distribution is dominated 
by the noise term in this region as seen in Fig. [1] (Left) , little information of the galaxy distribution from this region 
can be expected. Therefore, we neglect the non-linear effect of the angular power spectrum of the galaxy distribution 
in this paper. 

Usually, the full Fisher matrix for joint experiment of galaxy survey and CMB is obtained simply by adding each 



Fisher matrices: K, 



pCMB 
ij 



pgalaxy 
ij 



This method, however, does not include all the available information because 
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it does not account for the cross-correlation of temperature- galaxy C; and lensing potential-galaxy C; . To use 
the angular power spectrum C; ^ and C; for the estimation, and make the most of the available information, we 
consider all of the conceivable cross-correlations and calculate the full covariance matrix, which in this case is 8 x 8 
matrix while we assume that C; and C; ^ are not correlated. Here, in order to account for the difference of the each 
survey area, and to investigate the significance of the cross-correlation signal, we calculate the full Fisher matrix of 
the following form. 
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FIG. 3; The signal-to-noise from C; \ C; ^ and C;'^*, for zo = 0.8 (thin) and 1.2 (thick), respectively, with Planck (top) and 
CMBPol (bottom). The survey area and total number of the galaxies are fixed to /sky ~ 0.10 and A'g — 10^, respectively. 

where "/skv'™'^ " represents the available sky fraction in the Fisher matrix with (Survey) = {galaxy, CMB}, and 
"C^^ , • • • " denotes the angular power spectra included in the Fisher matrix with X,Y = {T, E, -0, g}. Case (i) does 
not include the cross-correlations between CMB and galaxy survey, C^ ^ and C; , and Case (ii) does not include the 
cross-correlations between temperature and galaxy or lensing potential, C^ ^ and C^ , while Case (iii) takes account 
of all cross-correlations. We shall compare the difference in these three cases in the following section. 

VII. RESULT 
A. Signal-to-Noise 

The signal-to-noise ratio (S/N) for a cross-correlation between X and Y can be estimated by |40l 



I) -fs.yf^i^l + 1)- 



{C'l 



{cr? + {Ci 



XX 



■Nr){cr 



■ND 



(40) 



We show the S/N value in Fig. |3] where we fixed the survey area as /sky = 0.10 and total number of the galaxies 
as A^g = 10^, respectively. The cross-correlations between temperature and lensing potential C; and temperature 
and galaxy Cj ^ are through the Integrated Sachs- Wolfe (ISW) effect imprinted in the CMB and the distribution of 
matter at late time. The ISW effect arise from the time- variation of the scalar metric perturbations and it is usually 
divided into an early ISW effect and a late ISW effect. The origin of the late ISW effect is from the time variation 
of the gravitational potential by the dark energy component and its effect emerges at low-multipoles. Therefore, the 
S/N both from Cj and Cj ^ saturate around /max — 40, although their amplitudes are different. 

On the other hand, the cross-correlation between lensing potential and galaxy has another feature. The survey 



which can explore the small scale region with large l-a 



can get large S/N from C; more than those from C; 



,T^|J 



and 



C; ^ while their amplitude is very small for the low resolution survey with small /max- The saturation of S/N from 
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C; at large /max for Planck in Fig. |3] (top two panels) is due to the noise contribution of the surveys. This justifies 
our omitting the proper modeling on small scales where non-linear evolutions are important. The signal-to- noise may 
be improved for high-precision future CMB survey, such as CMBPol, because it will be obtain much more information 
about small scale region than Planck. Since the cross-correlation between lensing potential and galaxy C; has larger 
S/N than other cross-correlation components, Cj would be more powerful tool when small-scale powers are observed. 
However, it should be noted that the correct non-linear model will be necessary on these scales. 

B. Parameter errors 

We show the la marginalized error of each parameter in Table ITVl and error contours in Fig. U) First, we compare 
the analysis methods, i.e., the Cases (i) - (iii) defined in Eqs. (|37| - (p9| . From the figure we find that the cosmological 
parameters, such as fl^, f^, and w are tightly constrained in Case (i), while the constraint on the non-Gaussianity 
parameter /nl is tighter in Case (ii) than Case (i). Because the main difference between Case (i) and (ii) is whether 
one includes Cj or C; , respectively, we can conclude that the cross-correlation between lensing potential and 
galaxy Cj is important to determine /nl rnore precisely. Assuming more high-precision CMB survey, CMBPol, the 
significance of C; for constraints on /nl are seen more clearly. 

Next, we investigate how the different galaxy sampling models affect the determination of /nl- For this purpose 
we consider two cases, the cases {a, (3) = (0.5, 3.0) and (a, /3) = (2.0, 1.5), and the results are shown in Fig. [5] The 
features of the two models are that the former has a gradual distribution and the latter has a sharp one around the 
peak redshift zq. Generally, the expected errors rapidly decrease with redshift zq for zq <\. As for the Case (iii) the 
error is almost independent of Zq for the case with (a, /3) = (2.0, 1.5), while it has strong dependence for the case 
with (a, (3) = (0.5, 3.0). However, translating the peak redshift zq to the mean redshift Zm, the tendency of the two 
models is similar to each other, although the constraint from the sharp distribution model is somewhat stronger than 
the gradual one, where the mean redshift z^ is determined by Eq. (|19p and the relation between z^ and zg depends 
on the galaxy sampling model. For the same value of the peak redshift zq, the gradual model represents relatively low 
mean redshift observation and the sharp model represents high redshift one. Therefore, the mean redshift of galaxies, 
rather than their distribution, determines how tight constraint one can obtain. 

The difference of the constraints due to sampling models can be attributed to the degree of the correlation between 
the lensing potential and the galaxy distribution. Because the sharp one has narrow peak around the peak redshift, 
it has much large correlation with the lensing potential in this narrow range. The degree of the correlation takes 
maximum value at certain redshift, and it gradually decreases above the redshift. This fact reflects that A/nl 
slowly increases with increasing of the peak redshift above the certain redshift. On the other hand, because the 
gradual distribution model has a wide peak, the galaxy distribution gas the lower degree of correlation with the 
lensing potential, even though the correlation exists over the wider range in fc-spacc. In other words, to constrain the 
primordial non-Gaussianity from galaxy-CMB lensing cross-correlation one should select the galaxies whose correlation 
with the CMB lensing potential becomes maximum. 

Finally, we focus on the observing redshift dependence for the constraints on the primordial non-Gaussianit /nl- We 
found that the constraints on /nl considerably depends on the mean redshift of the observations z,n, which is related 
with model parameter zq by Eq. (J19p . We show its redshift dependence in Fig. [51 The errors of /nl rapidly decrease 
with redshift at low redshift, zq < 1, while it increases slowly at high redshift, zq > 1- The effect of the primordial 
non-Gaussianity through the scale-dependent bias becomes large at high redshift, so that this is the reason why the 
smaller error of /nl can be obtained when the higher redshift is probed. On the other hand, the auto-correlation 
signal of the galaxy distribution becomes small with increasing redshift. This is an opposite effect to that from the 
scale-dependent bias for constraint on /nl and this is the reason why the error of /nl gradually increases at higher 
redshift, in particular, in Case (i) There are two reasons for increasing of the error of /nl at high-redshift. One is 
due to galaxy sampling model defined by Eq. ([5]). This analytic form may drop the information of the low-redshift 
galaxies at high zq. The other is that the cross-correlation between lensing potential and galaxy becomes weak at 
high redshift zq, as seen in Fig. [2] for /nl = (solid line). 

Moreover, we compare the constraints from various galaxy survey conditions and linear bias parameters fog in Fig.[Sl 
From this figure, we clearly see that both "depth" and "width" for galaxy survey arc important for constraint on /nl 
and what should be stressed is that the case of large bias b^ constrains more strictly than the case of low bias. The 
objects with large bias are affected by the primordial non-Gaussianity more strongly than the objects with low bias, 
so one of the key points to put constraint on the primordial non-Gaussianity is to explore the highly biased objects. 
This results indicate that some galaxy survey exploring the highly biased objects could constrain /nl in less than 10 
even with Planck, for example, A/nl ~ 5 at z^ = 2.0. 

The estimations given above do not take into account conditions of realistic observations, because we vary only the 
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FIG. 4: \a confidence limits on the pair (/nl, Q-i) in our 13-diniensionaI model. We show for Case (i) (dot-dashed line), Case 
(ii) (dashed line) and Case (iii) (solid line), respectively. We assume Planck (thick line) and CMBPol (thin line) for CMB 
survey, and for the galaxy survey the sky coverage and total number of galaxy sampling of /sky =0.1 and A'g = 10®. The 
parameters of galaxy sampling model are fixed as a = 0.5, /3 = 3.0 and zo = 1.8. 



peak redshift zo fixing the survey area /sky and the available galaxy samples A'g. In fact, in the real galaxy survey the 
survey area and the available galaxy samples will also vary due to the change of the observed mean redshift because 
the total observation time is finite as explained in Sec. IIVI Accounting the realistic galaxy survey condition, what 
strategy should we develop for constraining the primordial non-Gaussianity, e.g. deep survey or wide survey ? In the 
next section, we search for the best condition for constraining /nl- 



C. For the Galaxy Survey with Modeling 



Assuming observations like HSC, we estimate a more realistic constraint on /nl using the survey model introduced 
in Sec lIVI We show the results in FiglT] and the relations between various survey parameters, namely, peak redshift 
^0, survey area /sky, number density of sampling galaxy til and exposure time texp with a fixed observation time in 
Table El 

In Fig. [71 we find that there is a minimal point for the error of /nl around Zm — 0.7 and in this point, the 
constraints on /nl are A/nl ^ 20 for 6o = 2 and A/nl ^ 10 for &o — 4- From this result, the target redshift 
we should observe is not deep enough, in which case the surface number density of the sample galaxies is relatively 
small, and the observation area can be wide, riL — 7.0 [arcmin"^] and /sky — 0.35 (at z,„ = 0.7). (see Table FVl ) For 
the question whether we should make wide or deep survey, the answer is that we should select the wide survey. The 
signal of the primordial non-Gaussianity in the scale-dependent bias 6(z, fc) is more significant in relatively large-scale 
regions, as seen in Fig. |2l The noise contributions in large-scale regions are dominant by the cosmic variance due to 
the finiteness of survey area. On the other hand, small scale region is dominated by shot noise related to the surface 
number density of the galaxy samples riL, although the signal of the primordial non-Gaussianity is not sensitive there. 
Therefore, we conclude that the better constraining the primordial non-Gaussianity, /nli prefers wide area survey 
to the deep survey. However, note that this conclusion is derived by assuming Planck and HSC experiments with 
a fixed observation time. Accounting for the redshift dependence of the effect of the primordial non-Gaussianity 
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FIG. 5: The comparison of the la error of /nl in the cases with and without cross-correlations. The lines with simbols are 
against peak redshift zo and thick dashed line corresponds to the Case (iii) against mean redshift Zm, respectively. The left 
panel is for the case a = 0.5, /3 — 3.0 and the right one is for the case a = 2.0, /3 — 1.5, in Eq.([5]). Planck is considered for the 
CMB experiment and the galaxy survey parameters are taken as /sky = 0.10 and A^g = 10^. 
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0.0237 


0.0237 


w 


0.651 


0.197 


0.191 


0.195 


0.195 


0.253 


0.269 


0.270 


0.187 


0.195 


0.195 


iVoff 


0.135 


0.116 


0.116 


0.116 


0.116 


0.116 


0.117 


0.117 


0.116 


0.116 


0.116 


a 


0.00841 


0.00755 


0.00749 


0.00752 


0.00754 


0.00761 


0.00767 


0.00769 


0.00745 


0.00751 


0.00753 


/nl 






105.7 


39.9 


27.0 


104.9 


38.3 


25.9 


100.8 


36.9 


25.1 


bo 






0.0682 


0.0718 


0.0931 


0.0703 


0.0826 


0.1076 


0.0662 


0.0707 


0.0911 



TABLE IV: Constraints on cosmological parameters from the Fisher matrix methods for each case. We show the marginahzed 
la errors for the eleven-parameter or thirteen-parameter models. We assume Planck for CMB survey. For galaxy survey we 



assume that sky coverage and total number of galaxy sampling are /sky 
model are taken as a = 0.5, /3 = 3.0 for all cases. 



0.1 and A'g — 10 . The parameters of galaxy sampling 



through scale-dependent bias, we should keep it in mind that the deep survey also becomes important to put a tighter 
constraint on /nl- 



VIII. SUMMARY AND DISCUSSION 



In this paper we have estimated the constraints on the cosmological parameters newly taking into account the cross- 
correlation between CMB lensing and galaxy angular distributions. In particular, we have focused on the constraint on 
the primordial non-Gaussianity through the scale-dependent bias b{z, k) and estimated how much the cross-correlation 
between CMB lensing and galaxy would improve the constraint by the Fisher matrix analysis. In order to make the 
most general Fisher matrix analysis with CMB and galaxy survey experiments, we have taken into account the all 
the auto- and cross-correlations available which is expressed by the 8x8 covariance matrix as Eq. p3|) . We have paid 
particular attention to the CMB and galaxy survey just coming up now, namely Planck satellite and HSC, and also 
to the future experiments like CMBPol, LSST and so on. Our estimations are mainly based on the Planck and HSC 
surveys, however we also show some cases for comparison when CMBPol or ambitious survey conditions arc assumed. 
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FIG. 6: The la error of /nl for wide galaxy survey with different mean redshift. We assume Planck for CMB survey, and 
galaxy sampling model parameters are adopted as a = 0.5, /3 = 3.0. The thin and thick lines correspond to &o = 2 and 6o = 4 
cases, respectively. 



As for the constraints on the conventional cosmological parameters, the improvement can not be expected very 
much from the simple estimate in which the Fisher matrices for the CMB and galaxy surveys are combined properly 
even if the cross-correlations are properly taken into account. However, focusing on the determination of ./nl, we 
found that the role of the cross-correlation between CMB and galaxy is important, especially the one between CMB 
lensing potential and galaxy C^ contributes the determination of /nl- 

We have estimated the constraints on the /nl for the coming experiments. First, we gave the rough estimation 
in cases where the survey area /sky and galaxy samples iVg are fixed and only peak redshift zq is varied. It was 
found that the keys for more strictly constraining the primordial non-Gaussianity arc observing higher redshift and 
larger biased objects. Second, considering the realistic observations with fixed observation time, we have estimated 
the constraint on /nl, adapting the galaxy survey model in Ref. [3j], which is scaled to the HSC survey. As a result 
we found on optimized target redshift to be zq ~ 1, which brings us to a conclusion that we should make a wide 
survey rather than a deep survey for constraining the primordial non-Gaussianity with HSC-like observation. This 
is because the effect of the primordial non-Gaussianity through the scale-dependent bias is significant on large scales 
rather than small scales. The large-scale region is dominated by the cosmic variance related to the survey area while 
the small-scale region is dominated by the shot noise term related to the number of the galaxy samples. Therefore, 
we should explore the wide survey area rather than observing a lot of galaxies. However, we should keep it in mind 
that high-redshift and highly-biased objects are much affected by primordial non-Gaussianity, so that the deep survey 
will be also essential, in the future. 

The constraints on the primordial non-Gaussianity expected from HSC survey with Planck are : A/nl ^ 20 for 
6o = 2 and A/nl ^ 10 for bo = 4. Slosar et al. [l4| obtained constraints on /nl for highly biased tracers using 
available luminous red galaxy (LRG) or quasar (QSO) data from Sloan Digital Sky Survey (SDSS) [4l[ and CMB 
data from WMAP 5, with the error of A/nl ~ 30. The reason why the combination of HSC and Planck observations 
does not make significant improvement over the current constraints is explained as below. We have considered only 
800 hours for HSC and normal galaxies. The SDSS data is obtained over longer period of time, and QSOs seem to 
be observed more than normal galaxies at high-redshift. These constraints are weaker than those expected with the 
CMB bispectrum constraints achievable with an ideal CMB experiment, A/nl ^ few (Ref. [42|, |43)- However, the 
constraint on /nl presented in this paper depends highly on the galaxy survey condition, i.e., survey area /sky, number 
density of sample galaxies til and observing peak redshift zq. It is found that with some galaxy survey with Planck 
/nl could achieve A/nl ^ 5 at Zm = 2.0 (Fig.IH]). Recently it is reported that an ambitious future galaxy survey (like 
the LSST survey); which provides large survey area of 30,000 deg^ and highly biased galaxy samples, can measure the 
primordial non-Gaussianity with the order A/nl ^ 2 — 5 [15| . The other method using the full covariance of cluster 
counts for Dark Energy Survey (DES) can yield A/nl ^1 — 5 [13 • In Ref. [l^ [l^ , they add the Planck and galaxy 
survey Fisher matrices, i.e., Fij = Ff^ ^^^ -\- F^^^ , and do not include the cross-correlation between CMB and galaxy 

survey, i.e., Cj ^ an Cj , so that more tight constraint on /nl may be expected with these cross-correlations. In any 
case, it is worth pursuing how well we can put a constraint on non-Gaussianity of the local-type from the large-scale 
structure because it contains information on non-Gaussianity at different epoch from CMB and thus the constraint 
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FIG. 7: The la error of /nl in the case with fixed sky coverage /sky =0.1 and total number of galaxy A'^g = 10^ (Left), and 
in the case with modeling for HSC survey (Right). For the CMB experiment we consider two cases. The thick lines show the 
errors against peak redshift zo and the thin hnes are against mean redshift Zm. 

through the scale-dependent bias will be an important cross check against the CMB bispectrum. 



zo 



A/nl 



Uh arcmm 



/e 



sky 



0.8 
0.9 
1.0 
1.1 
1.2 
1.3 
1.4 
1.5 



0.51 
0.58 
0.64 
0.70 
0.77 
0.83 
0.90 
0.96 



30.1 
24.5 
20.7 
18.5 
18.6 
21.9 
29.1 
40.5 



0.9 
1.9 
3.7 
7.0 
12.4 
20.9 
34.0 
53.5 



0.42 
0.41 
0.40 
0.35 
0.25 
0.14 
0.06 
0.02 



TABLE V: The relation of parameters with fixed observation time, between the peak redshift zo, the mean redshift Zn 
la error of /nl, the number density of sampling galaxy riL and the survey area /sky, for the HSC-like survey {bo — 2). 
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